Power corrections in e+e tt+tt , K K and B — > Ktt. tttt 
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CLEO-c measurements of the timelike form factors F^, Fk at ^/s — 3.671 GeV provide a direct 
probe of power corrections (PC's) at energies near ms. PC's in -FV./c and B Kit , tttt are separated 
into perturbative parts, and soft parts obtained from fits to data. In the latter are 0(10) 

larger. Soft dominance of same order is expected in the QCD penguin PC's, and is required by the 
B Kn,mr fits. The CP asymmetries ^^^^o, Cn^^o are well determined, agree with experiment, 
and provide a clean test for new physics. Greater PC soft-overlaps for tttt, and much larger ones for 
PP than VP, VV shed light on the special nature of the pseudo-Nambu-Goldstone bosons. 
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Much effort has gone into the theoretical description of 
B decays into light meson pairs. Apart from being of in- 
terest in QCD, the issue has important implications for 
new physics search strategies which rely on comparing 
decay rates and CP asymmetries in different final states. 
The decay amplitudes can be organized into expansions 
in powers of 1/mh. The leading power (LP) contributions 
are calculable in QCD factorization (QCDF) [ij in terms 
of universal non-perturbative quantities. Numerous lead- 
ing power predictions for B M1M2 decays are in gross 
conflict with the data. In i? ^ Kn, tttt the direct CP 
asymmetry A^+^- is too small, Aj^+^- is too small and 
of wrong sign, Ai^+^- « Aj^+^a , contrary to observation, 
and the branching ratios Br^o^o, Br^o^o are too small. 
A possible explanation is that certain PC's are of same 
order as or larger than their LP counterparts and have 
large strong phases, due to non-perturbative effects. 

Continuum e+e^ — > AI1M2 light meson cross sections 
at y/s 3.7 and 10.58 GeV at the charm and B factories 
provide a direct probe of PC's in the timelike vector- 
current matrix elements {MiM2\q"f^q\0) ■ Perturbative 
calculations of PC's on the light-cone contain IR log- 
divergent terms of the form as(/ih)(l/\/s)" ln™(Y^/A), 
signaling the breakdown of short /long-distance factoriza- 
tion. A represents a physical IR cutoff on the longitudinal 
momentum of, e.g., a valence quark in the convolution in- 
tegrals of light meson light-cone distribution amplitudes 
(LCDA's) with hard scattering amplitudes. We therefore 
divide the PC's into perturbative and non-perturbative 
parts (soft overlaps), where the former are defined by im- 
posing A, fifi > 1 GeV, and fits to the data give ranges 
for the latter. For example, the vector-current form fac- 
tor PC's are written as SF = SF^"'^- + (5F" p-. In this 
letter we focus on PC's in e+e^ — > tt+tt^, K^K^ and 
B Ktt, tttt. Details, and results for VP, VV final 
states will be given elsewhere 

The continuum timelike form factors Fk , at ^/i = 
3.671 GeV are il 



0.075 ±0.009, I = 0.063 ±0.004. 



(1) 



The calculable LP contributions arise at twist-2 in the 
LCDA's, and fall like 1/s [J. We obtain 
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FIG. 1: 5F^^^' vs. A. Solid curves are for Qs and central val- 
ues of LCDA parameters evaluated at = A. Inner bands 
(blue) are for variations of LCDA parameters. Outer bands 
(yellow) include variation of € [Max[l, A./2\,^/s\. Er- 
rors are added in quadrature. Dashed lines outline the 'outer 
bands' for asymptotic LCDA's. 



where the errors due to variation of the first two LCDA 
Gegenbauer coefficients Q, and the scale at which they 
and as are evaluated, y/s to -s/s/S, are added in quadra- 
ture. Thus, the form factors must be dominated by PC's. 
A dual-resonance model calculation oi F^^ |6|, which in- 
cludes the first three p resonances explicitly, is consistent 
with the twist-2 result for ^/s > 3 GeV 17|. Evidently, 
duality has set in, implying sensitivity to PC's. 

5Fk,t! enter at or twist-4 perturbatively, and to 

first approximation fall like 1/s^. We obtain SF^"^' from 
convolutions of two twist-3 valence quark LCDA's with 
the tree-level hard-scattering amplitudes (twist-4 valence 
quark LCDA's contribute negligibly). The model param- 
eters of [5| are employed for the LCDA's. Perturbative 
higher Fock state effects are of same twist and order of 
magnitude, and therefore would not alter our result for 
the soft PC's. Fig. 1 shows the ranges obtained as the 
cutoff A in the divergent terms is varied from ^/s to 1 
GeV. A is roughly the lowest gluon virtuality allowed. 
Ranges for asymptotic LCDA's are also shown. There 
are large accidental cancelations between asymptotic and 



non-asymptotic effects at lower A. The magnitudes of 
each separately therefore give a better indication of the 
size of perturbative effects for A near 1 GeV. Compar- 
ing the asymptotic plots to Eq. [1] it is clear that the 
dominance of the PC's in Fk,tt is due to their soft parts, 

^i^n p /^i^pcrt. ^ 0(^0), SF'^P /SPP"''- = 0(10). (3) 

Similar soft enhancement would account for Fjr{mj/^) « 
0.10, as extracted from J/^" decays [3]. 

At LP, the form factors obey canonical SU{3)p fla- 
vor symmetry breaking, i.e., {F^/Fk)lp ~ f^/fx ~ 
0.7. However, |exp. = 1.20 ±0.17, implying that 

\6F^/SFk\ > 1, e.g., (1.9 ± 0.36) for construc- 
tive interference between LP and PC effects. Apparently, 
the soft-overlap is greater for tttt than KK, providing a 
source for large SU{3)f breaking. 

The SU{3)p diagrammatic representation gives a con- 
venient general classification of the B — > Kt:, tttt ampli- 
tudes 0. The B Ktt amplitudes are 
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The CKM factor Ap = VpbV*^, and there is a sum over p — 
u,c. The B- {B°) Br's are given by |Ap {tbo /tb-\A\^) 
in our normalization. T [ai) and C (02) are the 
color-allowed and color-suppressed 'tree' amplitudes. 
PP (04,6), Pew («7,9), and P^^ (as, 10) are the QCD, elec- 
troweak penguin (EWP), and color-suppressed EWP am- 
plitudes, respectively. In tttt the corresponding quantities 
are primed ('). These amplitudes consist of LP parts 
(the corresponding QCDF coefficients ai [l| are given in 
parenthesis) Tlp, Clp, etc., and PC's 5T, 5C, etc.. Pew 
and A are EWP and 'tree' weak-annihilation (WA) PC's. 

The Kit PC fit is dominated by {SPPle^^", \6C\e'^'^ , 
\ST\e^^'^ . The strong phases are defined relative to the 
corresponding naive-factorization amplitudes. Allowing 
for a soft contribution that is an order of magnitude 
greater than its perturbative part leaves A negligible. 
The electromagnetic u, c-loop penguin contractions are 
negligible in the EWP PC's, so we drop their 'p' super- 
scripts. The approximate relation (up to an unknown 
SUi3)F breaking factor) SP^y^ « 3/2C9/C1 5C, and the 
order of magnitude relations (equalities in the one-gluon 
exchange approximation) SPew ^ 3/2Cio/CiSC and 
SP§^ ~ 3/2C9/C2A relegate the EWP PC's to a mi- 
nor role. The Ci are Wilson coefficients for the AB = 1 
effective Hamiltonian [Tj . In terms of the measured un- 
certainties, the impact of (5Pew on ^^k+tt- j Br^+j^o , 




FIG. 2: Diagrams for perturbative B Mi A/2 PC's. Crosses 
indicate all places where the gluon can attach. The Qi are 
the AB — 1 effective Hamiltonian operators 

Ak+tt- is found to be ^ 2.5a, 1.5a, 1.5a, respectively. 
The impact of the remaining EWP PC's and A on all 
observables is < la. The tttt PC fit is dominated by 
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The 'tree' WA ampli- 



tude E plays a significant role in tt+tt , and the tttt EWP 
PC's are negligible. 

AU Br's and Ax+tt-, Ak°tt+i Ak+-k°i ^tt+tt-j 'S'tt+tt- 
are re quir ed to lie within their la HFAG experimental 
errors ;iOl|. The LP amphtudes are evaluated in QCDF 
to NLO 1] . The NNLO corrections llj would not have a 
substantial impact on our fit results. The LP inputs are 
varied uniformly within their errors. The Wilson coeffi- 
cients, as, and the LCDA parameters [5| are evaluated 
at the scale /i^ G [to?, /2 , m^] . The fit allows 5P^ and 
5P^ (JP'" and SP"^) to differ substantially, as suggested 
by their perturbative contributions, see Fig. 4c. In fact, 
|jp«/^pc| < 2, 3 for 69%, 87% (ISP"" /SP'^l < 2,3,4 for 
55%, 81%, 92%) of the Kn (tttt) scatter points in Fig. 3. 
According to Figs. 3a, c \5P'^\ ^ \PIp\, and SC dominates 
in C with a significant strong phase 6c ~ 30°. The fit 
results for C/T in Fig. 3d allow for smaller magnitudes 
than obtained in SU{3)f fits, see e.g. [ij]. The need for 
a large strong phase difference Sc — St is well known. The 
SM predictions for the time-dependent CP asymmetries 
Sk^^o e [0.62,0.86] and C^.^o = -A^^^G [0.07,0.18] 
are consistent with the measured ranges [10|, see Fig. 3e. 
Nevertheless, the measured uncertainties exceed those for 
the fits, making these observables a good place to look 
for new physics, see e.g., [l^. Note that the Kn results 
show little dependence on \6P'^ /5P'^\. 

The TTTT P''^ fit in Fig. 3b is multiplied by fK/f-n for 
comparison with Fig. 3a. This time, |<5P'^| Z |PlpI- 
Canonical SU{3)f breaking at LP gives /k/ZttPlp ~ 
Pip. However, fK/f-n^P"^ appears to be larger than 
5P'^, in accord with SF^j > (5_FV.This conclusion is rein- 
forced if \5P'^ /5P"^\ is bounded from above, e.g. < 3, 
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FIG. 3: K-K, TTTT fits: (a) SP" (blue), P£p (red), P" (green), (b) 
SP'" (yellow) and for |5P'"/5P"=i < 3 (blue), P{3, (red), P'" 
(green), (c) 5C (blue), Clp (red), C (green). Strong phases 
are relative to naive factorization (d) |C/r| (blue), \C' /T'\ 
(green) vs. their strong phase differences (e) SM predictions 
for Sji^^o, Cj^^^o and HFAG averages. 



see Fig. 3b. C /T' is plotted in Fig. 3d. The lower bound 
on its magnitude increases from roughly 0.15 to 0.25 for 
|^P'"/(5P'^| < 3. The predicted ranges for A^+^o, A^o^o 
are [-0.06, +0.06], [-0.95,0.55], respectively, show little 
dependence on \5P'^ /5P'^\, and are consistent with the 
HFAG averages 0.06 ± 0.05, 0.43 ± 0.25. 

The perturbative B — + Kit, tttt PC's are obtained from 
the diagrams of Fig. 2. They depend on two renormal- 
ization scales: (i) linked to the energy release of the 
decay, as in the LP amplitudes, and (ii) /x/j, linked to the 
IR cutoff A. The Wilson coefficients and one as factor 
in Figs. 2b, d (associated with the w,c -loops) are eval- 



uated at ^f). The other as factors and the LCDA pa- 
rameters are evaluated at jih- The Wilson coefficients in 
Figs. 2a,c,e are NLO, and Ci is LO in Figs. 2b,d. We have 
checked that the loop diagrams in Fig. 2 eliminate the 
dominant leading log^b scale dependence (c>c Cias/T^) 
in ^pp.poit-. Products of twist-2,3 x twist-3 K,tt valence 
quark LCDA's are included in the convolutions with the 
hard-scattering amplitudes. Again, higher Fock-state ef- 
fects would not alter our conclusions. 

Our results are summarized in Fig. 4. The largest real 
and imaginary contributions to ,5P'''^^'P®'^*' come from 
Fig. 2b. The QCD dipole operator Q^g (Fig. 2e) and 
QCD penguin WA [Q^^..^ in Fig. 2a) contributions are a 
factor of 2-3 smaller in magnitude, and real. All three are 
dominated by diagrams in which a gluon is not attached 
to the B. They contribute at twist-3 (l/m;,). (Other 
contributions to (5p(')p ^poit.^ g_g^ Q3,..,6 in Fig- 2c and 
Fig. 2d, are negligible.) In the limit A Aqcd they 
would give rise to two energetic outgoing quarks, and 
soft-overlaps for both light mesons, as indicated by a 
quadratic dependence on log ms/A. This is also the case 
for SF^°^', which depend quadratically on log ^/s/A. 
Therefore, the relative importance of non-perturbative 
PC's in JP(')P should be similar to Eq.El given the prox- 
imity of energy scales. Indeed, comparison of Figs. 4a,b,c 
with the fit results in Figs. 3a, b implies 

gpin c. n.p./^p(/) c, pert. ^ 0(^Q) 

(5(7(') pert, jg (j^g j-fj pjg^ 2c. The diagrams again 

enter at twist-3 (1/mf,). In the A — > Aqcd limit only 
the spectator quark light meson would be produced via 
a soft-overlap, as indicated by a linear dependence on 
log ms/A. Therefore, we expect ,5C(') ■ " p/JC^') ■ to 
be smaller than the QCD penguin ratio. According to 
Figs. 3c, 4d (for low A) this may well be the case, but 
the errors are too large to reach a definitive conclusion. 

We emphasize that the similarity of Eq. p to Eq. [3] 
supports the importance of the charm-loop |14| . dipole, 
and WA PC's in the full QCD penguin amplitudes, and 
increases our confidence that the B — > Kir, tttt 'puzzles' 
can be accounted for as in the fits. Given the penguin 
PC soft dominance, the origin of fa/ fir SP''^ > SP'^ is the 
same as SF^^ > SFk, a larger tttt soft-overlap. That cer- 
tain soft PC's are of same order as their LP counterparts 
should be viewed as accidental, rather than a breakdown 
of the 1/E expansion. The latter would be signaled by a 
violation of power counting among the PC's themselves. 

The asymptotic LCDA plots in Figs. 4a, b imply that 
the charm quark loops can yield large strong phases in 
the full penguin PC's. Their presence may be indi- 
cated by the fits for in Figs. 3a, b. However, a siz- 
able strong phase could also originate from soft-overlaps 
in light meson production, independently of such loops, 
explaining the large strong phase in C. The existence 
of such a mechanism would be confirmed by a non- 
vanishing strong phase difference between the LL and LT 
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FIG. 4: B — * Ktv, tttt perturbative PC's vs. A. Solid curves 
are for /ij, = mi,, fih = A, central values for other inputs. In- 
ner bands (blue) are for fih = A, fJ.b & [mb/2 jTUb], other 
inputs varied within uncertainties. Outer bands (yellow) in- 
clude variation of fj.h G [Max[l, A/2], ms]. Errors added in 
quadrature. Dashed curves outline 'outer bands' for asymp- 
totic LCDA's. Inner brown, outer orange bands in (a) are for 
Qsg- Black (brown) curves in (c) are for /k / /tv SP"^^^K 

(i=longitudinal, T^transverse) polarization amplitudes 
in e'^e~ — > p'^p~ at the T(4S'). This can be searched for 
with the current B factory data sets 0, ■ 

The asymptotic LCDA plots for ^p'^^.P^rt. ^^d 
gp'c, pert, jj-^ -pjg^ support the possibility, allowed by 
the PC fits, that the total up and charm penguin PC's 
can have a large hierarchy, e.g., (5P^'^ "/(JP*-'-* ;6 3. The 
perturbative analysis also implies that they can have 
large strong phase differences. Similar considerations in 
B — > (jjK imply that a large negative shift in the time- 
dependent CP-asymmetry relative to sin 2/3 is possible 
e.g., « -(0.10 -0.15) i. 

The QCD dipole operator contribution to JP^^p^''* is 
also plotted in Fig. 4a. Eq. O then implies that the total 
Qsg PC could be of order Plp, see Fig. 3a. It is notewor- 
thy that this is an order of magnitude larger than its LP 
matrix element. Therefore, the size of new CP violating 



dipole operator contributions currently probed could be 
0(10) smaller than expected based on LP analyses. 

Continuum e+e^ — > VP, VV cross sections at ^/s « 3.7 
and 10.58 GeV and penguin-dominated B VP, VV 
data can be accounted for if soft-to-perturbative PC ra- 
tios vary from 0(1) to O(few) rather than 0(10) as 
for PP. This appears to shed light on the "p — tt puzzle" 
(and the smaller "K* — K puzzle"), i.e., how to simulta- 
neously understand the tt as a state like the p, and a 
nearly massless Nambu-Goldstone boson. In a physical 
picture of the tt (and A') that addresses this puzzle [3| , 
the valence q and q are always a distance < 1/Aqcd 
apart. However, r^r ~ r-p « 1/Aqcd is accounted for 
by a much larger soft cloud of higher Fock states, at- 
tributed to a zitterbewegung-like motion of the tightly 
bound valence quarks. This may account for the large 
PP soft-overlaps, as well as the larger tttt soft-overlaps 
because the tt approximates a Nambu-Goldstone boson 
more closely than the K. 

Note added: After completion of the work reported 
here we were informed by L. Silvestrini that his group 
has also carried out a Ktt PC fit, with similar predic- 
tions for SksTt°, Cfc^Tr", which has since appeared IITI . 
For an earlier discussion by them, see M. Pierini in ll_4|. 
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